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Abstract 

In a remarkable numerical analysis of the spectrum of states for 
a spherically symmetric black hole in loop quantum gravity, Corichi, 
Diaz-Polo and Fernandez-Borja found that the entropy of the black 
hole horizon increases in what resembles discrete steps as a function 
of area. In the present article we reformulate the combinatorial prob- 
lem of counting horizon states in terms of paths through a certain 
space. This formulation sheds some light on the origins of this step- 
like behavior of the entropy. In particular, using a few extra assump- 
tions we arrive at a formula that reproduces the observed step-length 
to a few tenths of a percent accuracy. 

However, in our reformulation the periodicity ultimately arises as 
a property of some complicated process, the properties of which, in 
turn, depend on the properties of the area spectrum in loop quantum 
gravity in a rather opaque way. Thus, in some sense, a deep explana- 
tion of the observed periodicity is still lacking. 

1 Introduction 



Recently, a large computerized analysis of the spectrum of states for a spher- 
ically symmetric black hole in loop quantum gravity was carried out by 
Corichi, Diaz-Polo and Fernandez-Borja [lUlHsHll. The analysis focused in 
particular on the entropy for the black hole. The theory l6ll7l[8ll9l[T0l[TTl[T2l 
predicts the area-dependence of the entropy as 
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y is the Barbero-Immirzi-parameter, and yc ^ numerical constanllJ. Fixing y 
to be equal to y c insures consistency with the Bekenstein-Hawking entropy 
law. The numerical analysis d 12113114] was confined to relatively modest 
black hole sizes, up to a few hundred Planck-lengths squared in terms of 
the black hole area. Still, it required the calculation of up to 1 0^^ states. 
The work beautifully confirmed the leading order behavior linear in A and 
it was even powerful enough to confirm the next to leading order contri- 
bution. But what is more, it discovered a completely unexpected periodic 
behavior of the entropy: There is a very distinctive periodic signal super- 
imposed onto the linear growth ([T]). Its period was found to be proportional 
to y, 

AA=yxli (2) 
where x was estimated in [2 1 to be 

X ~ XCDF = 8.80. (3) 

This signal leads to a stair-case like behavior for the entropy as a function 
of the area and was hence called entropy quantization in ||T| . 

The observation of this period is remarkable in several ways. For one, ^ 
can be read as a sort of effective equidistant quantization of area [IJ. Thus 
although the area spectrum in loop quantum gravity is not equidistant but 
rather consists of sums of terms of the form 

Aj = STryl^yjOTTT, i G N/2. (4) 

^ does provide a point of contact with the ideas of Bekenstein (see for ex- 
ample millll) on area quantization. This point gets even more pronounced 
when one takes into account that, as first observed in fT], the value Xcdf is 
very close to 81n(3) ^ &.7?>&?>9?>. Based on a heuristic quantization of the 
black hole and the Bekenstein-Hawking-relation for entropy and area, an 
equidistant area spectrum had been conjectured ||14| , with 

AA = 41n(k)l^ (5) 

and k an integer. Later, k = 3 was suggested by Hod fTS] using a corre- 
spondence to the frequencies of quasinormal modes of the black hole in a 

Vc depends on the precise definition of which states are counted as surface states. There 
are two possibilities, leading to two slightly different values for yc : Just counting the m- 
and b- labels (ex. (8l|9|) leads to Yc ~ 0.237, whereas also counting the j-labels (ex. |12|) 
leads to Yc ~ 0.274. 
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certain limit. Dreyer IIT6II observed that in loop quantum gravity one could 
reinterpret ^ as the smallest nonzero area eigenvalue, making a certain 
choice for the parameter y that seemed natural at that time, but it was later 
realized that that y would not lead to the Bekenstein-Hawking-relation for 
entropy and area. Given all this, it is intriguing that an equidistant spec- 
trum with A A involving ln(3) seems to reappear here. We should not fail to 
mention however, that there is an additional factor of 2y in ((2]) as compared 
to (O, so ultimately it is not clear wether there is a deep correspondence 
with the ideas of Bekenstein and Hod. For more comments on this see HJ, 
for some considerations of physical consequences we refer to lH. 

Secondly the phenomenon seems to be robust against the way the the split 
between surface and bulk degrees of freedom is done. As indicated in 
AA just depends on y, and x (and not on y^), and is hence independent of 
the way the states were counted]! 

Finally, there does not seem to be a straightforward way to explain this phe- 
nomenon from the theory. After all, as (H) shows, the area spectrum is quite 
complicated, non-equidistant, and does not, in any obvious way, determine 
the constant X- We should note that the Aj of (|4]) do become approximately 
equidistant for large j (for details see [ISJ), so it might at first seem that this 
gives an explanation for the observed periodicity. However, the spacing in 
(HJ becomes close to multiples of 47rylp « 12.56 yip. Comparing this with 
(lU, it is clear that this can not furnish a simple explanation for the value of 
X- 

With the present paper we aim to contribute to an explanation of this phe- 
nomenon of entropy quantization. We use two main ideas. The first is to 
reformulate the combinatorial problem of enumerating the physical states 
for the black hole horizon in terms of paths built from a set of elementary 
steps. The second idea is to use a statistical description of the set of paths, 
very similar to a random walk. 

These ideas, together with some assumptions on the numerical distribu- 
tion of steps will enable us to calculate an approximation to the period AA 
that reproduces the observed step-length to a few tenths of a percent accu- 
racy. AA arises as some sort of 'resonance' in the area spectrum (H)). While 
we think that this is a nice result, it is not a complete analysis and expla- 

^It is only a posteriori that one might decide to adjust the freely specifiable parameter y 
to take the value yc and hence make {l) consistent with the Bekenstein-Hawking entropy 
formula. 



3 



nation of the phenomenon. To start with, since we have to make some 
assumption and approximations, our result for AA is not exact, and its un- 
certainty is hard to determine. Thus we are unable to confirm or rule out 
that X = 8 In (3). Furthermore our analysis will not determine wether the 
phenomenon will persist for arbitrarily large areas. We will discuss these 
points in more detail, below. 

After the present work was finished, we became aware of results ||T7| by 
Ansari. His is a very nice analysis of the spectrum of the full area opera- 
tor, breaking it down into an infinite set of equidistant sub-spectra. For the 
black hole horizon, only a subset of that spectrum is relevant, due to the 
horizon boundary conditions, and the exclusion of edges that are a sub- 
manifold of the horizon, as well edges that pierce the horizon from within 
the black hole. Ansari's methods break down for this "reduced" area oper- 
ator. Still it is not inconceivable that his results are connected to the entropy 
quantization, and thus may represent another way of looking at the prob- 
lem tackled in the present work. 

Finally we would like to bring to the reader's attention that there is in- 
dependent very interesting work on the way [5J which, using somewhat 
similar methods, will shed more light on the phenomenon of entropy quan- 
tization. 

The paper is organized as follows: In the next section we review the com- 
binatorial problem of enumerating horizon states of the black hole and for- 
mulate it in terms of steps and paths. Section|3]is concerned with statistical 
considerations and with a computation of X- We finish with a discussion of 
the results and future perspectives in SectionHl 

2 Counting states by counting paths 

In the present article, we will not review any of the physics behind the 
description of an isolated horizon in loop quantum gravity. We refer the 
interested reader to IZllHl- Here it will suffice to spell out the combinatorial 
problem to which finding physical states of the black hole horizon is ulti- 
mately reduced. The quantity of greatest interest to us is N (I), the number 
of horizon states that are eigenstates of horizon area with eigenvalue in the 
interval I. As we have said in the introduction, there are two ways to count 
states, depending on where one draws the line between bulk and bound- 
ary degrees of freedom. Both lead to qualitatively similar results for N(I). 
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Here we will consider only one of these ways, the one that was laid out 
in IHl in detail. In |9| this problem was revisited and given a very simple 
formulation. It was shown that N(I) is the number of ordered sequences 
(rai)t, TRi G Z*/2 such that 

^ mi = and STtyl^ Y_ \/|TTLiKimi| + 1) G I. (6) 

i i 

The cormection to the quantum geometry of the horizon is that sequences 
(rai)i with these properties are labels of physical states of the horizon. Let 
us simplify even further and get rid of all the units, by defining n(a) = 
N (Syrylpa). It will also be useful to introduce the shorthand 



a(m) = V|m|(|m| + 1), m G Z/2. 

Finally we take from I^.TO] the idea that the counting problem can be sim- 
plified by implementing the two conditions of (H)) in separate steps. We 
define 



[m^,m2, . . .), mi G Z*/2 : ^ mi = j, ^ a( 



mi 



It will be instructive to reinterpret the state labels as paths. To this end, 
introduce the space § = R+ x Z/2. Let us call a sequence of points (pi)i in 
this space a path, and the differences pi+i — pi the steps of the path. Now let 
us call a path (pi)i allowed if 

• it starts at 0, i.e. pi = (0, 0), and 

• the steps pi_|_i — piare of theformv(mi] = (a(mi],mi) for some mi G 
Z/2. 

Obviously, then, we can associate to a state labeled by (mi , . . . vnn) the al- 
lowed path (0, v(mi ), v(mi ) + v(m2), . . .). 

What we are really interested in is the number n(I) of states with the area 
in an interval I. In the language of paths this is given as 

n(I) = the number of allowed path with endpoint in I x {0}. 

Given any number R > there is a finite number of allowed paths that end 
within [0, R] X Z/2 C S. It is easy to enumerate them, and we have written 
a little Mathematica routine that does this for us. In Figured] we show the 
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Figure 1: Paths that end below R = 4 (left) and below R = 8 (right). 



results, for R = 4 and R = 8, by plotting, for each path, its endpoint [a,j) as 
a dot in the diagram. In effect these diagrams contain all information about 
the functions n(I) and n(a, j). The most obvious feature of the results is the 
striking regularity that they exhibit. We should mention that many of the 
points plotted in Figured] lie on top of each other. For R = 8 there is for ex- 
ample a total of 76619 points to plot. So the regularity in the concentration 
of points becomes even more striking if one plots the density of endpoints 
of paths. We have done so in Figure |2] for R = 8. 

Now an important question is the following: Is the period we see in these 
figures the one that was observed in fTJ|2H3l|4l? The answer is yes: The a = 
const, lines drawn in the Figures [T] and |2] are regularly spaced at intervals 

Aa = 0.35. 

Their correspondence to the period exhibited in the data is clear. More- 
over, physical states are associated to the points on the line M+ x {0} in the 
diagrams. Thus the periodicity on that line corresponds directly to the pe- 
riodicity seen in IIl|2l|3l|4j. We have also plotted density on this line and 
logarithm of density in a sliding interval [a — Aa/2, a -I- Aa/2] (Figure |3l 
The latter starts to show the staircase of IllOISjIlI near a = 8. 

It must be said that the amount of data that we have assembled is much 
smaller than that handled in the much more sophisticated analysis lU |2l 
Based on our data alone it would be premature to conclude that a 
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Figure 2: Density (left) and logarithm of density (right) of paths that end 
below R = 8 (in fiducial units). 



ln(density+ 1 ) ln(density+ 1 ) 




Figure 3: Logarithm of density of physical states, in a small interval around 
a (left) and in an interval adapted to the periodicity (right). 
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periodicity in the entropy is present. We are however confident that our 
plots show the onset of the pattern that [11121 EHH has demonstrated much 
more clearly and to much higher values of area. 

What we have done so far certainly does not amount to an explanation 
of the periodicity. We merely reformulated the problem of enumerating 
surface states of the black hole into one concerning paths in the space §. 
We then observed that the periodicity found in (ij [H |3l HI does apparently 
not just govern paths corresponding to physical states (the line j = in the 
figures) but a larger class of paths. In the next section we will show that 
using the image of steps and paths can be very helpful in the analysis of 
the pattern. 

3 AA and the statistics of the steps 

Now that we have exhibited the pattern in our reformulation through paths 
and steps, let us return to its explanation. The reformulation can shed 
new light on the issue as follows: Imagine for a moment that all the steps 
that would be allowed in allowed paths were just integer multiples of one 
basic step. Then regularity in a diagram like Figure [T] would obviously 
result. This is not the case for the system at hand: The allowed steps 
v(ra) = (a(ra), ra) are not integer multiples of one another. Moreover, even 
if the steps were approximately multiples of one basic step (as could be ar- 
gued is the case for the area spectrum |[T8|), stringing together many steps 
would in general lead arbitrarily far away from points in a regular pattern. 
So at first sight the consideration of a situation with just integer multiples 
of a basic step does not seem to lead anywhere in terms of explaining the 
pattern observed here and in lTl|2l|3l|4J. We should however keep in mind 
that what we want to explain is not a completely rigid phenomenon. It is 
not so that there are no states that fall outside the pattern observed in Figure 
[H it is rather that the majority clusters around some evenly spaced points. 
Moreover the observed pattern is really the result of thousands (or in the 
case of |[ll|2l [31111, more like of 1 0'^'^) points, which in our language are each 
obtained by taking many steps. So what we should be looking for is not a 
formula that describes all the details of the spectrum, but something that 
explains why it is statistically likely for a point to lie in one of the clusters. 

Coming back to the steps and paths, our idea is as follows: If we can 
demonstrate that the steps are multiples of a single step on average in a suit- 
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able sense and moreover that the variance of the steps around this average 
is very small, then we can at least explain that a pattern formed at low areas 
will reproduce itself for some time. Let us try to make this argument more 
precise. 

Let us assume that there is a well defined probability distribution p(m) for 
the occurrence of a step v(m) in a path corresponding to a physical state. 
Let us furthermore assume that we can treat the individual steps in a given 
path corresponding to a physical state as independently distributed and 
with the distribution p(ra), to a good approximation. Then let us write 

a(m) = I(m)Aa + e(m) (7) 

where I(ra) shall be an integer and e(m) < Aa. Now we consider a path 
with n steps v(rai) that starts somewhere on the lattice AaZ x Z/2. The 
distance of the endpoint of this path to the corresponding lattice point is 

n 

6(n,{v(mi)}) = Y_ e(TrLi). 

i=l 

Now we look at this quantity under the probability distribution. Because 
of our assumptions we can use the central limit theorem to approximate 

(5(n)) « n(e(m)), {8[nf - {8[n)f) « n{e(mf - {e[m)f), (8) 

where the averages on the left of these approximate equalities are expec- 
tation values in the ensemble of physical paths with n steps, whereas the 
averages on the right are in the ensemble of steps, 

(f)= n^nUm) 

TaGZ*/2 

for f a function on Z*/2. (|8]) shows that if we can choose Aa in lO such that 
(e(m)) = then we can expect the path to remain near the lattice as long as 



n(e(ra)^) < Aa, or n< 



[Aa 



,2 



Let us apply this reasoning to the problem at hand. What we need is (a) 
information about the statistics of the steps involved, and (b) we will have 
to make an ansatz for the function I(ra) of dZ]). 

As for (a), luckily we have at least some information on the statistics of the 
horizon states. In [9J it was shown that upon picking a path at random out 
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Figure 4: The possible steps starting at (left) and a set of these steps su- 
perimposed onto the spectrum of paths (right) 

of all physical paths ending below some ao, the probability to find v(ra) as 
the first step is 

p(m] f« exp (^-27ryM\/!m|(lm| + 1)^ 

with ~ 0.2375 the numerical constant that makes the sum of the p(ra) 
over ra equal 1, and the approximation good as long as a(m) is small com- 
pared with ao- What we would rather like to know is a different probability, 
namely that of finding v(m) as first step in a path among all the paths ofn 
steps and ending below ao- Now, as long as n is large, (but not as large as 
ao), the dependence of this probability on n should be rather weak, and 
thus we assume that it is proportional to the p(ra) above. Thus we will 
work with an ensemble of steps with the above probability distribution, 
i.e. we will define the average for a function f on Z*/2 by 

(f)= Y. f(m) exp (-27rYMv/|m|(|m| + 1)) . 

rn.GZ./2 

Now we turn to part (b): We have to make an Ansatz for the function I of 
(IZ|. To that end, we inspect Figure S] which shows the basic steps, and how 
they fit into the pattern of the allowed paths. One sees very clearly that 
the regularity in the pattern of the paths and of the steps are related. More 
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precisely a(m+ 1 /2) — a(m) « 3Aa/2. Moreover there is a shift of one unit, 
independent of ra. Altogether, we will write 

a(ra) = ^^-2111+1^ Aa + e(ra). (9) 

which in turn defines the quantities e(vn.), once Aa is fixed. 

Now we can proceed as outlined above. We want to determine Aa, and we 
want to do it in such a way that ( e ( ra) ) is zero, otherwise any pattern would 
be washed out. Taking averages of (O and using the condition (e(ra)) = 
indeed determines Aa: 

3(ra) + 1 

This in turn fixes 

(a(m)) 



e(ra) = a(ra) — [3vn + 1 ' 



'3(m) + T 

Numerical evaluation of these formula can be done very easily on a com- 
puter. We find 

Aa« 0.34952, (e(ra)^) w 0.00019156 (10) 

Let us put these numbers in perspective and into context. First of all, we 
find that the standard deviation for the e(ra) is very small: 

25 



That means that only after a number n of steps of the order of 625 (= 25^) 
do we expect to deviate from the pattern substantially, as we have argued 
before. This means that at the very least our results are significant for the 
black holes of small area as considered here and in (TH2H3l|4|. Secondly, 
our value for Aa compares nicely with the one obtained in HHlHSllll. Let 
us illustrate this in terms of the parameter x- It is related to our Aa by 
X = SttAu, which compares to the result xcdf of [1. 2, 3» 4J as follows: 

X~ 8.7843, XCDF -8.80 XCDP-X ^ o_oo^29, 

Xcdf 

so we agree with the quoted reference to within a fraction of a percent. 
What is more, we seem to be even closer to the conjectured value 8 ln(3): 

81n(3] « 8.7889, ^^T^i^)-X _ 0.00053. 

X 

We will discuss these findings further in the next section. 
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Figure 5: The set of paths ending below R = 8, with a distorted area spec- 
trum a'(m) 

4 Discussion and outlook 

What we have done in the present paper is to give an explanation of the 
phenomenon of entropy quantization in loop quantum gravity by means 
of a formulation using paths, steps, and their statistics. The periodicity in 
the spectrum of horizon states arises as some sort of resonance (O in the area 
spectrumll That ^ works so well has to do with the fact that the area spec- 
trum is nearly equidistant, v^(j-M) ~ j + 1 /2. The precise relation Q that 
gives small (e(ra)^) (and hence the value of A A) does however depend on 
the details of the area spectrum as well as on the quantum boundary condi- 
tions (i.e. the j = constraint). Thus it is intimately related to properties of 
the area quantization of loop quantum gravity, however in a rather opaque 
way. A nice way to confirm this is to redo the calculation of allowed path, 
however using a slightly distorted area spectrum a'(ra) = a(ra) + 1 / (I Ora). 
One can see that the details of the result (Figure |5]l change quite drastically 
as compared to the undistorted spectrum. One does however still recog- 
nize a lot of regularity in the result. Our interpretation is that because the 
area spectrum is still approximately equidistant, one again sees regularities 

^We should stress that although we talk about a resonance, and our model for the black 
hole involves some sort of random walk, these are not physical processes. They have nothing 
to do with the dynamics of the black hole. 
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emerge, whereas the precise pattern has changed because the resonance 
condition to achieve small (e(ra)^) is now different. 

Our explanation seems to be quite successful quantitatively, as we recover 
the results of lU 121 |3l HI for x- However here already one problem of our 
approach becomes apparent: Several approximations go into the determi- 
nation of our value for X/ and we have little idea how accurate the result 
actually is. 

As for other aspects of the phenomenon, some can be explained by our ap- 
proach, while others remain mysterious. In particular, we want to remark 
the following: 

(1) The phenomenon occurs for both ways to count [H HI El IH but here we 
have considered only one. We note however that the other way of counting 
states (i.e. the inclusion of j-labels) can be seen as a refinement of the count- 
ing we have done here. To be more precise, each path corresponding to a 
physical state that we have counted here corresponds to one or more phys- 
ical states as counted in the other scheme. Thus the pattern that we have 
observed is bound to appear also in the other scheme, possibly modulated 
further by some other effects that come from the details of the counting of 
the j-labels. Thus it seems to us that the explanation of the phenomenon 
given here also applies to the other counting scheme. 

(2) It was observed fH |2l |3l IH that the phenomenon goes away when not 
implementing the condition that ^|rai = 0. From our Figures [T] and IH 
as well as our expression for the area spectrum (HJ it appears that the pat- 
tern is shifted by A a/2 between lines with 1] even and lines with 2j odd. 
This explains at least why the pattern gets washed out considerably when 
summed over all ] . It could be, however, that a (substantially weaker) pat- 
tern with a spacing of Aa/2 remains, and it would be interesting to look for 
it in numerical data. 

(3) It was observed [Tl|2l|3l|4l that x is very close to 81n(3). The situation 
here is very tantalizing, in that on the one hand our result for x moves even 
closer to the conjectured value. On the other hand, since our treatment is 
only approximate, we can not draw any conclusion from this. 

(4) It has been conjectured [H 121 |3l lU that the phenomenon continues to 
be present even for macroscopic black holes. From our treatment it does 
seem that the pattern should start to get washed out once the black hole 
is so large that the dominant paths are longer than about 600 steps. We 
do not understand the mechanism at work in the generation of the paths 
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well enough to present this as a result, however. Rather, we must leave this 
question open for future research. 

Altogether, we think that the approach taken affords interesting insights, 
but it does not give answers to some crucial questions. Moreover it may 
even be questioned wether the present approach can answer these ques- 
tions even when worked out in more detail, since it uses statistics and 
some heuristics. Therefor it would be interesting to pursue alternative ap- 
proaches. One possibility that comes to our mind is to analyze in detail the 
properties of the Laplace-Fourier-transform of n(Q, j) as determined to a 
good approximation in [10] (and which can be determined exactly, as far as 
we can see, with similar methods). Its structure, and in particular its poles, 
should contain information on periodic phenomena of n(a, j) along lines 
j = const. We will pursue this approach elsewhere. 
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